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Abstract. Results for the n + N ^ A,!^ + K reactions in nuclear matter of 
Rcf. [0 are presented. To evaluate the in-medium modification of the reaction 
amplitude as a function of the baryonic density we introduce relativistic, mean- 
field potentials for the initial, final and intermediate mesonic and baryonic states 
in the resonance model. These vector and scalar potentials were calculated using 
the quark meson coupling model. Contrary to earlier work which has not allowed 
for the change of the cross section in medium, we find that the data for kaon 
production at SIS energies are consistent with a repulsive i<'+-nucleus potential. 



1. Introduction 

The properties of kaons in nuclear matter have recently attracted enormous interest 
because of their capacity to signal chiral symmetry restoration or give information 
on the possibility of kaon condensation in neutron stars |^]. Studies with a variety 
of models [|[ ^ indicate that the antikaon potential is attractive while the kaons 
feel a repulsive potential in nuclear matter. However, the analysis of available 
data on production from heavy ion collisions at SIS energies ||, || contradicts 
the predictions that the kaon potential is repulsive. The comparison between the 
heavy ion calculations and the data [sj indicates that the if^-meson spectra 

are best described by neglecting any in-medium modification of the kaon properties. 
Furthermore, the introduction of even a weakly repulsive K^-micleus potential results 
in a substantial underestimate of the experimental data on kaon production. 

Since in heavy ion collisions at SIS energies ||^, ^ the X^-mcsons are 
predominantly produced by secondary pions, we investigate the kaon production 
reactions, n + N —^ A, E -I- K , in nuclear matter. We calculate the in-medium reaction 
amplitudes, taking into account the scalar and vector potentials for incident, final and 
intermediate mesons and baryons. 



2. Mean-field potentials for mesons and baryons 

We use the quark- meson coupling (QMC) model which has been successfully 
applied to the various problems of nuclear physics [LlJ and the studies of meson and 
hyperon properties in a nuclear medium [l^ . The Dirac equations for the quarks 
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and antiquarks in the hadron bags {q = u,u,d or d, hereafter), neglecting the Coulomb 
force, are given by Q: 

[ij ■ dx - ms]ips{x) (or ^s{x)) =0. (2) 

The mean-field potentials for a bag in nuclear matter are defined by — g^a, — 
g^LO and = g'^b, with g^, g^ and the corresponding quark-meson coupling 
constants. 

The hadron masses in symmetric nuclear matter relevant for the present study 
are calculated by: 

_ (n, -I- nq)ni + (n, + nj)^?, - Zh ^ i^j^*3^ ^^^^ 9ml 



Rl ' 3 ' dRh 



= 0, (3) 

where fi* — ^ x'^+{Rlm*y , with m* = mq—g^a and — \/xf+iR^jns)^ , and 

Xq and Xs being the bag eigenfrequencies of the corresponding quarks. In Eq. (^), Uq 
(riq) and Ug (n^) are the lowest mode light quark (antiquark) and strange (antistrange) 
quark numbers in the hadron, h, respectively, and the Zh parametrize the sum of the 
center-of-mass and gluon fluctuation effects, and are assumed to be independent of 
density. The parameters are determined in free space to reproduce their physical 
masses. We chose the values — 5 MeV and rrig = 250 MeV for the current quark 
masses, and Rn = 0.8 fm for the bag radius of the nucleon in free space [Q. We 
stress that only three coupling constants, g^, g^ and g^, are adjusted to fit nuclear 
data - namely the saturation energy and density of symmetric nuclear matter and the 
bulk symmetry energy. Exactly the same coupling constants, g%, g^ and g^, are used 
for the light quarks in the mesons and hyperons as in the nucleon. However, for the 
kaon system it was phenomenologically necessary to increase the strength of the vector 
coupling to the non-strange quarks in the , and thus, we will use the stronger vector 
potential, lA^VJ, for the if+-meson [Q. The scalar {U!}) and vector (L/^') potentials 
felt by the hadrons, h, in nuclear matter are self-consistently calculated by: 

U': ^Us^ ml - m,„ = {nq - nq)V^ - hV^ , (FJ ^ 1.4 for K+), (4) 

where, /s is the third component of isospin projection of the hadron, h, and the p 
meson mean field potential, V^*, is zero in symmetric nuclear matter. Then, within 
the approximation that the mean field potentials are independent of momentum, the 
four-momentum of the hadron is modified by, = [yP^+rn^ + U^,p), which 
modifies not only the kinematical factors such as the flux, the phase space and the 
threshold, but also modifies the reaction amplitudes. 

Furthermore, we include the in-medium modification of the resonance masses, 
which appear in the reaction amplitudes in the resonance model. In view of its 
numerous successful applications elsewhere, we base our estimate on the QMC 
model [0, n2i, and use those estimated values IpJl for the in- medium resonance masses. 



3. TT + N A + K in nuclear matter 

Now we apply the resonance model Q to calculate the in-medium amplitudes 
focusing oin:-\-N—fA + K. We consider kaon and hyperon production processes in 
ttN colhsions shown in Fig. 0. (See Refs. H, || for tt -I- ^ S -I- i^T.) 
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Figure 1. A" and A production processes included in the resonance model 



We extend the model by including medium modification of the hadron properties, 
not only in the kinematic factors such as the flux and the phase space, but also in the 
reaction amplitudes, as discussed in the previous section. (See also Ref. 

We found that the maximal downward shift of the reaction threshold in nuclear 
matter occurs at baryon densities around pB — O.Gpo (po = 0.15 fm~^), and the 
maximum of the downward shift of the tt + A^^ A, E + if reaction threshold amounts 
to roughly 30 MeV. We also found that at baryon densities pB > 0.2 fm~'^ thresholds 
for the TT + N ^ Aj'E + K reactions are higher than those in free space. 

We show the energy dependence of the total tt~ + p ^ A + cross section in 
Fig. ||. The calculations for free space are in reasonable agreement with the data, 
as shown by the solid line. The dashed line in Fig. || shows the results obtained for 
nuclear matter at pb = Po, while the dotted line is the calculation at pB = 2po- 
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Figure 2. Energy dependence of the total cross section, n~ + p — > A + K^, as 
a function of the invariant collision energy, y^, calculated for different baryon 
densities. The data in free space are taken from Ref. mm. The lines indicate our 
results for free space (solid) and for nuclear matter at baryon density pg = po 
(dashed) and ps = 2po (dotted) (po = 0.15 fm"^). 



4. Impact on heavy ion studies 

It is expected that in relativistic heavy ion collisions at SIS energies nuclear matter 
can be compressed up to baryonic densities of order pB — 3po Hi- The calculation of 
the time and spatial dependence of the baryon density distribution is a vital aspect 
of dynamical heavy ion simulations. The baryon density pB available in heavy ion 
collisions evolves with the interaction time, t, and is given by the dynamics of the 
heavy ion collision. In the following estimates we investigate the density dependence 
of the production cross section for central central heavy ion collisions. 

To calculate the X+-meson production cross section averaged over the available 
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density distribution we adopt the density profile function obtained by dynamical 
simulations ||l5| for Au + Au collisions at 2 AGeV and at impact parameter = 0. 
This can be parametrized as 

PB(t) = p-max exp {[t - t ]VAt2) , (5) 

where the parameters, Pmax — 3 po, f = 13 fm and At = 6.7 fm, were fitted to the 
heavy ion calculations ]15[. 




(GeV) 

Figure 3. Energy dependence of the total cross section, ■k~ + p ^ h + , as 
a function of tlie invariant collision energy, -y^. The data in free space are taken 
from Ref. |14[ . The solid line indicates our calculation for free space. The dashed 
line shows the cross section calculated by averaging over the density function 
profile lli given by the time evolution obtained for Au + Au collisions at 2 AGeV 
(sec Eq^l)). 

The total cross section for the tt^ + p ^ A + K'^ reaction integrated over the 
time range 5 < i < 23 fm and weighted by the time dependent density profile given 
in Eq. (||), is shown by the dashed line in Fig. ^. The limits of the t integration were 
taken from the simulations of the Au + Au collision time evolution in Ref. |^ . The 
circles and solid line in Fig. || show the experimental data in free space and the 
calculations in free space, respectively. 

One can see that the total cross section averaged over the collision time (time 
dependent density profile) for the + p A + reaction is quite close to the 
result given in free space integrated up to -^s ~ 1.7 GeV. Indeed, the results shown 
in Fig. ^ actually explain why the heavy ion calculations with the free space kaon 
production cross section might quite reasonably reproduce the experimental data. 

Furthermore, the total cross section averaged over the time dependent density 
profile, shown by the dashed line in Fig. ||, should additionally be averaged over the 
invariant collision energy distribution available in heavy ion reactions. The number of 
meson-baryon collisions, NmB, for the central Au + Au collisions at 2 AGeV is given 
in Ref. [ p6[ as a function of the invariant collision energy, ^/s. It can be parametrized 
for ^/s >1 GeV as 

dN,nBld^s = No exp {[^fs ~ ^ ]V[AV^]'), (6) 

where the normalization factor — 6x10^ GeV^^, while ^/s^ — 1 GeV and A^/s = 
0.63 GeV. Note that, at SIS energies NmB is almost entirely given by the pion-nucleon 
interactions, and heavy meson and baryon collisions contribute only to the high energy 
tail of the distribution in Eq. (^) - with quite small densities 16 . 
average the calculated, in-medium, total cross section for tt" + p 
by the dashed line in Fig. 0, over the available energy distribution given in Eq. 



Finally, if we also 
A + if" , shown 
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we obtain an average total kaon production cross section of <K> — 65 /ib for central 
Au + Au collisions at 2 AGeV. This result is indeed compatible with the calculations 
using the free space total cross section of the 7r~ +p — > K + K^ reaction, which provide 
an average total kaon production cross section of <K> = 71 /ib for central Au + Au 
collisions at 2 AGeV. Note that the inclusion of even a slight modification of the K'^ 
mass because of the nuclear medium (without the corresponding changes introduced 
here) leads to a substantial reduction of the inclusive spectra (by as much as a 
factor of 2 or 3), compared to that calculated using the free space properties for the 
relevant hadrons 

We stress that at SIS energies reaction channels with a S-hyperon in the final state 
play a minor role, because of the upper limit of the energy available in the collisions. 
The downwardly shifted tt + N ^ + K reaction threshold at small densities is still 
quite high. 



5. Conclusion 



Our present study shows that if one accounts for the in-medium modification of the 
production amplitude (i.e., the in- medium properties of the i^T^-meson and hadrons) 
correctly, it is possible to understand production data in heavy ion collisions at SIS 
energies, even if the i<r+-meson feels the theoretically expected, repulsive mean field 
potential. The apparent failure to explain the production data if one includes the 
purely kinematic effects of the in- medium modification of the ii'^-meson and hadrons, 
appears to be a consequence of the omission of these effects on the reaction amplitudes. 
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